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Abstract 

Nonleptonic B to charmonium decays generally show deviations from the 
factorization predictions. For example, the mode B~ — ► K~XcO has been 
experimentally observed with sizeable branching fraction while its factorized 
amplitude vanishes. We investigate the role of rescattering effects mediated 
by intermediate charmed meson production in this class of decay modes, and 
consider B~ — ► K~h c with h c the J PC = 1^ cc meson. Using an effective 
lagrangian describing interactions of pairs of heavy-light Qq mesons with a 
quarkonium state, we relate this mode to the analogous mode with XcO m the 
final state. We find B(B~ — * K~h c ) large enough to be measured at the B 
factories, so that this decay mode could be used to study the poorly known 
h r .. 



PACS number: 13.25.Hw 
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I. INTRODUCTION 



The precise test of the Standard Model description of CP violation in the B sector 
is among the most challenging efforts pursued at present experimental facilities. It goes 
through the measurement of many observables, such as CP asymmetries and B meson 
branching fractions which are sensitive to CKM angles. In order to extract meaningful 
information from experimental data, a reduced theoretical uncertainty is required, and this 
is a particularly demanding task in the case of nonleptonic B decays for which a completely 
reliable and general computational scheme has still to be developed. 

For two-body nonleptonic B decays, which concern us in the present paper, the determi- 
nation of the transition amplitude reduces to the calculation of the following matrix element 
of the effective Hamiltonian governing B — > MiM 2 [1]: 



In (1.1) \ are CKM matrix elements, q(//) Wilson coefficients evaluated at the scale /i and 
Oi a set of four-quark operators. So, neglecting corrections to the r.h.s. of eq.(l.l) that 
are suppressed by inverse powers of Mw, the analysis of the decay amplitude involves the 
calculation of hadronic matrix elements of four-quark operators. The oldest prescription, 
which could be used to evaluate any generic form (1.1), is the naive factorization ansatz 
that expresses the matrix elements of four-quark operators as products of hadronic matrix 
elements of quark currents. 

Let us consider B~ — > K~ M £c which is pertinent to our discussion; M Sc is a meson 
belonging to the charmonium system. Neglecting the annihilation term which is suppressed 
by the CKM factor V u b, the effective Hamiltonian Hw driving the decay reads as 



h w = -^{v^cMOM + cMom) - VttV£5>M<W} + h - c - ( L2 ) 



A(B -> M 1 M 2 ) 



G p 
V2 



Y^Kc^M.M^O^B) 



(1.1) 



where 



(cb) V -A(sc)v- A 



o 2 



(sb) V - A (cc) V -A 




^4(6) 



(sibj)v-A^2(qjqi)v-A[v+A] 



(1.3) 



^7(9) 



■z(sb)v-AY, e q {qq)v+A[v-A] 




-(Sibj) V -A e q(qjqi)v+A[V-A] 
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(i and j are color indices and (qq)v=pA = ?7 M (1 TTs)?)- The corresponding expression of the 
factorized amplitude is: 

„4 F (£T - K-M Sc ) = ^y cb T/*[a 2 (//) + ]T a i (/i)](X-|(sfo)y_ A | J B-)(M ac |(c C )y^|0) 

i=3,5,7,9 



(1.4) 



where are combinations of Wilson coefficients: a 2 = c 2 + and Oj = q + with N c 
the number of colors. 

Eq.(1.4) shows the drawbacks of the naive factorization approach: first, the scale and 
scheme dependence of the Wilson coefficients Cj(/i) is no longer compensated by a corre- 
sponding dependence of the hadronic matrix element, and secondly, the product of hadronic 
matrix elements does not contain any strong phase. 

Great amount of work has been done since this formulation of factorization has been 
put forward, aiming either at finding alternative procedures or at changing the ansatz itself. 
An improvement consists in adopting a generalized factorization ansatz, with the Wilson 
coefficients eij(/x) replaced by effective (process independent) parameters a| ff to be fixed using 
experimental data. In some cases this method reproduces the correct order of magnitude of 
the branching ratios [2]. Other methods, such as QCD-improved factorization [3], PQCD 
[4], SCET [5], QCD sum rules [6,7], can only be applied to selected classes of nonleptonic 
transitions. 

In B to charmonium decays, generalized factorization indicates the existence of sizeable 
nonfactorizable contributions. For example, the experimental branching fraction BiB — > 
K~ J/ip) can be fitted using jaf^ = 0.2 — 0.4 depending on the B — > K transition form 
factor which parametrizes the matrix element (K~\(sb)v-A\B~) in (1.4) 1 ; |a| ff | = 0.38 ± 
0.05 is obtained using the form factor in [8]. This must be compared to the value a 2 = 
0.163(0.126) computed for fn^rrib) = 4.4 GeV and = 290 MeV in the naive dimensional 
regularization (or 't Hooft-Veltman) scheme [1], a value which does not change significantly 
by varying m b (m b ) and A|^. The difference between af 1 and a 2 witnesses the presence of 
nonfactorizable effects in this decay mode. 

However, the most compelling evidence of deviation from factorization comes from the 
observation of B~ — > K~Xco, with Xco the lightest cc scalar meson. The measured branching 
fraction is: 

B(B~ - K- Xc0 ) = (6.0^-1 ± i.i) x 10- 4 (1.5) 
B(B~ -> K- X co) = (2.4 ± 0.7) x 10~ 4 (1.6) 

for BELLE [9] and BABAR [10] Collaborations, respectively. While the experimen- 
tal amplitude evidently is non- vanishing, the factorized amplitude (1.4) is zero because 
(Xco\(cc) VtA \0) = 0. Interestingly, the decay occurs at a rate comparable to B~ — > K~J/ip 



Since the other Wilson coefficients are numerically small, one can safely consider only the con- 
tribution proportional to a 2 . 
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since, for example, ggL^g-Jfcjg = (0.60±g;f| ± 0.05 ± 0.08) as reported by BELLE Collab- 
oration [9]. 

Analyses of the two modes B~ — > K~XcO, K~ J/ip in the framework of QCD-improved 
factorization show that perturbative QCD corrections are not able to reproduce the experi- 
mental branching ratios, giving either small contributions or producing infrared divergences, 
a signal of uncontrolled nonperturbative effects [1 1] . 

In ref. [12] we investigated the possibility that the deviation from the factorization pre- 
dictions in B — > charmonium processes may be ascribed to rescattering processes, essentially 
due to intermediate charm meson exchanges represented by diagrams of the type depicted 
in fig. 1. Rescattering effects in heavy meson decays have been considered recently, for ex- 
ample, to explain the observation of some OZI-suppressed decays of -^(3770) [13], or as 
possible contributions to B -> tttt [14], B -> K^ir [15], [16], B s -> 77 [17]. We found that 
rescattering effects could be sizeable, enough to produce a large branching ratio as observed 
in B~ -> K'Xco- 

Here we wish to reconsider the problem, since other decays modes have vanishing fac- 
torized amplitude [18] and can be used to test the rescattering picture. One of them, 
B~ — > K~h c with h c the lowest lying J PC = l + ~ cc state, deserves particular attention. 
The meson h c was searched [19] and observed [20] in pp annihilation, and searched in p — Li 
interactions [21]; the reported mass and widths are m^ c ~ 3526 MeV and T < 1.1 MeV. 
It is listed by the Particle Data Group among the particles requiring confirmation [22]. If 
B~ — > K~h c proceeds with a sizeable rate, this decay could be used to study the properties 
of h c by looking either at its hadronic transitions: h c — > J/ipn , p°ir°, /i 1 / (980), h^KK, 
. . . , or at its radiative decay modes: h c — > 7^7, Xcol, etc. 

This paper is devoted to such an investigation. Moreover, it aims at improving the 
analysis of rescattering effects in B to charmonium transitions reducing the dependence of 
the rescattering amplitude on unknown hadronic parameters, such as the strong couplings 
among different mesons. We introduce an effective lagrangian describing the interaction of 
all the low-lying £ = 1 charmonium states to pairs of open charm mesons, based on the 
spin symmetry for the heavy quark in the infinite heavy quark mass limit. This allows to 
express all the couplings in terms of a single hadronic parameter, as shown in Section III. 
A similar relation is derived for the couplings of £ = cc mesons to pairs of D^). Using 
such relations it is possible to analyze various rescattering amplitudes; their calculation is 
reported in Sections II and IV, while the conclusions concerning the possibility of using B 
decays to study the h c are drawn in the last Section. 

II. MODEL FOR CHARMED MESON RESCATTERING CONTRIBUTIONS 

As for B~ — > K~XcOi the factorized amplitude A F (B~ — > K~h c ) in (1.4) vanishes since 
the matrix element {h c \{cc)v^A\^) is zero due to conservation of parity and charge conjuga- 
tion. This does not imply that the decay is forbidden, as other decay mechanisms can be 
invoked, namely h c production via cc pair creation in the color octet configuration. From 
the hadronic point of view, one can also consider the decay as proceeding by rescattering 
processes induced by the same (bc)(cs) effective weak Hamiltonian in (1.2), processes that 
essentially account for a rearrangement of the quarks in the final state. Such effects are not 
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CKM suppressed, and their role must be assessed by explicit (even though model depen- 
dent) calculations. Notice that color octet and rescattering descriptions can represent two 
ways to describe the same physics underlying the nonleptonic transition, looking from the 
short- distance or the long-distance view points, respectively. 

We consider rescattering processes corresponding to the decay chain B~ — > X? C Y" 5 ~ — > 
K~M 5c , where X and Y are open charm resonances primarily produced in weak B~ tran- 
sitions. The lowest lying intermediate states and Y" 5 ~ are the mesons £)(*)- and £>W°, 

(*) 

and we describe their rescattering by the exchange of D)( resonances, as depicted in fig. 1. 

In order to analyze the diagrams in fig. 1 we need the weak vertices B — > D^D^ and 
two strong vertices, one describing the coupling of a pair of charmed mesons to kaon, the 
other one representing the interaction of the charmonium state h c to a pair of mesons. 
All nonperturbative quantities entering in such vertices can be related to few hadronic 
parameters once the infinite heavy quark mass limit is adopted. 

In the following Section we analyze the couplings of the charmonium states to pairs of 
open charm mesons. Here we consider strong interactions of mesons Hq containing a single 
heavy quark Q which can be described in the framework of the Heavy Quark Effective Theory 
(HQET) [23], exploiting the heavy quark spin and flavour symmetries holding in QCD for 
itlq — > oo. In this limit the heavy quark four velocity v coincides with that of the hadron 
and it is conserved by strong interactions [24]. Because of the invariance under rotations 
of the heavy quark spin sq, states differing only for the orientation of sq are degenerate 
in mass and form a doublet. When the orbital angular momentum of the light degrees of 
freedom relative to Q is £ = 0, the two states in the doublet have spin-parity J p = (0~, 1~) 
and correspond to (-D( s ), D*^), (B^ : B*^). This doublet can be represented by a 4 x 4 
matrix: 

Ha = (^f) " M a75 ] , (2.1) 

with corresponding to the vector state and M to the pseudoscalar one (a is a light 
flavour index). The fields M a and M* contain a factor ^Jm M (*), with m the meson mass. 

In the infinite heavy quark mass limit it is possible to express weak as well as strong 
matrix elements involving heavy mesons in terms of few universal quantities. Let us con- 
sider the weak amplitude B~ — > for which there is empirical evidence that the 
calculation by factorization reproduces the main experimental findings [25]. Neglecting the 
contribution of the operators O3-10 in (1.3) we can write: 

(dv-dW\h w \b-) = ^v cb v; sai (D ( * )0 \(v - Ay\B-)(DV-\(V - A),\0) (2.2) 

with ai = Ci + C2/N C . In the infinite heavy quark mass limit, the matrix elements in (2.2) can 
be written in terms of a single form factor, the Isgur-Wise function £, and a single leptonic 
constant F [23]. The B~ — > matrix elements read: 



< D°{v')\V^\B~(v) > = y/m B m D £{vv'){v + v'y 

< D*°(v',e)\V^\B-(v) > = -i^m B m D * £(v ■ v') e* p e^v^ (2.3) 

< D*%v\e)\A»\B~{v) > = y/m B m D * £{v ■ v') e*J(l + v ■ v')g^ - lA/"] , 
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v and v' being B~ and _DM° four-velocities, respectively, e the D* polarization vector and 
£(v ■ v') the Isgur-Wise form factor. The weak current for the transition from a heavy to 
a light quark Q — > q a , given at the quark level by <? a 7 M (l — 7s) Q, can be written in terms 
of a heavy meson and light pseudoscalars. The octet of the light pseudoscalar mesons is 

iM 

represented by £ = e / , with 



M 



7i- - v^ + v h k° 



(2.4) 



2 

and / ~ /„- — 131 MeV, and the effective heavy-to-light current, written at the lowest order 
in the light meson derivatives, reads: 

K = ^Tr[7"(l - 75)^U • (2-5) 

In this way the matrix elements < 0|g a 7 M (l — ^^D^iy) >, defined as 

< 0\q a Yl5c\D a (v) > = f Da m D y 

< 0\q a Yc\D* a (v, e) > = f D *m D ^ (2.6) 

can be related to the single quantity F since fr> a = fo* = 

It is also possible to write down an expression for the strong couplings involving heavy 
mesons and the kaon. The D^D^K couplings, in the soft px — > limit, can be related to 
a single low energy parameter g, as it turns out considering the effective QCD lagrangian 
describing the strong interactions between the heavy D^D^ mesons and the octet of the 
light pseudoscalar mesons [26]: 

d = ig Tr[H bl ^AlH a ] (2.7) 

with the operator A given by 

A l * a = l(t%t-td ll $) ha (2.8) 

and H a = 7°#]7°. This allows to relate the D^D^K couplings, defined through the 
matrix elements 

<D (p)K-( q )\D*-(p + q : 6)> =g DrD o K - (e-q) 
<D*°(p,r ] )K-(q)\D:~(p + q,e)> = i e a ^ p a epq^l g D .- D « K - , (2.9) 



to the coupling g: 



g 

9d*-d°k- = 2 v /m D m D * — 

JK 



9d* s -d*°k- ~ - 2 \ — —4- ■ (2-10) 

V m D* JK 

All the above expressions are valid in the infinite limit for the charm quark mass. We neglect 
corrections due to the finite mass of the charm quark. 
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III. COUPLINGS OF PAIRS OF HEAVY-LIGHT MESONS TO QUARKONIUM 

STATES 




The other strong vertex in the diagrams in fig.l involves h c and a pair of open charm 
mesons. Also in this case we exploit the infinite heavy quark mass limit. For mesons 
with two heavy quarks QiQ 2 heavy quark flavour symmetry does not hold any longer, but 
degeneracy is expected under rotations of the two heavy quark spins. This allows us to build 
up heavy meson multiplets for each value of the relative angular momentum t. For I = 
one has a doublet comprehensive of a pseudoscalar and a vector state, r\ c and J/ip in case 
of charmonium. The corresponding 4x4 matrix reads as [27]: 

(Y)^-^(V)- < 3j > 

with L M = J/xjj and L = r] c in case of cc. For i = 1, four states can be built which are 
degenerate in the heavy quark limit. The corresponding spin multiplet reads: 

where, in the case of cc, X2 = Xc2, Xi — Xci an d Xo — Xco correspond the spin triplet, while 
the spin singlet is hi = h c [28]. Also the fields in (3.1), (3.2) contain a factor y/m, with m 
the meson mass. 

Using (3.1) and (3.2), together with (2.1) representing the heavy-light Q\q a pseudoscalar 
and vector states, it is possible to write down the expressions for the effective couplings 
between heavy-heavy mesons and pairs of heavy-light mesons we are interested in. For 
i = 1 Q1Q2 state, the most general lagrangian describing the coupling to two heavy-light 
mesons Qiq a and q a Qi can be written as follows: 

& = i^-Tr [P^^H 2a {n lllx + Q 2 v^)H la ] + h.c. + (Q 1 <- Q 2 ) (3.3) 

where Oi and are two coefficients, H\ a is given in (2.1) and H 2a is the matrix describing 
the heavy-light mesons with quark content q a Q2- 

H 2a = [M>^ - M: 75 ] {^f^ . (3.4) 

Due to the property P^v^ = only the term proportional to Qi contributes, and therefore: 

A = i^-Tr [P (Q ^H 2al ,H la ] + h.c. + (Q 1 » Q 2 ) , (3.5) 

where g± — g\ • Qi. This expression accounts for the fact that the two heavy-light mesons 
are coupled to the heavy-heavy state in S-wave, and therefore the matrix elements do not 
depend on their relative momentum. Moreover, this expression is invariant under indepen- 
dent rotations of the spin of the heavy quarks, representing the decoupling of the spin in the 
infinite heavy quark mass limit. This can be easily seen considering that under independent 
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heavy quark spin rotations Si G SU(2)q 1 and S 2 G SU(2)q 2 the following transformation 
properties hold for the various multiplets: 



H\ a — > SiH la Hi a — > Hi a S\ 
H 2a —> H 2a S 2 H 2a — > S 2 H 2a 

piQiQz) _^ $ iR (QiQ2) ftiQiQ*) R (QiQ2) gt 

Eq.(3.5) shows that a unique coupling describes the P^HH interaction, i.e. the same cou- 
pling controls the interaction of heavy-light mesons both with the three Xc states, both with 
h c . In particular, from (3.5) it follows that: 

< £( 8 )(pi,ei)I?( s )(p2)|/ic(p,e) > = gD* s) D (s) hM ■ e ) 
< D *s)(Pi, £i)D* s) (p 2 , e 2 )\h c (p, e) > = i g D * s) D* s) h c ^raP a ^ f e* T ~e* 2 a (3.7) 

with 



9D* s) D (3) h c = -2gi^rn hc m D(s) m D * s 



rn 



D 



(3.8) 



Analogously: 



< J D( : s )(Pi ) ei) J D( s )(P2,e 2 )|Xco(p) > = -9D* (s) D* (s)Xc0 ( e *i ' e 2) 

< D W (pi)D W (p2)|Xco(p) > = -to (s) D (3)Xc0 (3.9) 



with 



2 

^?.) D (V<" = -^9iy/™^rn D . a) 

9d m d (s)Xc0 = -2V3g lv ftf^m D{s) . (3.10) 

The subscripts (1) and (2) refer to the meson with a charm and an anticharm quark, respec- 
tively; e, ei and e 2 are polarization vectors. 

Eqs.(3.7)-(3.9) show that spin symmetry produces stringent relations between the cou- 
plings of Xco an d h c to open charm mesons, relations that we exploit below. Moreover, 
they also imply that the couplings of a single charmonium state to open charm pseudoscalar 
and vector mesons are related in absolute value and in sign as well, a property that al- 
lows a proper analysis of the amplitudes in fig.l where the relative signs between different 
amplitudes play an important role. 
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For the I = states represented by the multiplet (3.1), the interactions with the heavy- 
light vector and pseudoscalar mesons proceed in P-wave and can be described by a lagrangian 
containing a derivative term: 



C 2 = ^Tr 
2 



R (QlQ2) H 2a H 



la 



+ h.c. + (Q 1 <-> Q 2 ) (3.11) 



which is also invariant under independent heavy quark spin rotations. The action of the 
derivative produces a factor of the residual momentum k, i.e. the quantity for which the 
hadron and the heavy quark four momentum differ: MhV^ = vtiqv^ + k^, k being finite in 
the heavy quark limit. The couplings of heavy-light charmed mesons to J/tp follow from 
(3.11): 

< Dl s) (p 1 ,e 1 )D* {s) (p 2 ,6 2 )\J/ip(p,e) > = 9d* s) d^ 

[(e-e*)(e* 1 . q )-(e. q )(e* 1 .e*) + (e.e* 1 M. q )\ 

< £( s) (pi,ei)D w (p2)|J/^(p,e) > = g D * 3)D(s)i ,ie^ a ye' I e* 1 a q T (3.12) 

< D {s) ( Pl )D {s) (p 2 )\J/iP(p,e) > = g D(g)D(g)i ,(e ■ q) 

where q is the difference in the residual momenta of the two heavy-light charmed mesons 
q — k± — k 2 . Since pi = m n (*)V + k\ and p 2 = m n (*)V + k 2 , then q = p 1 — p 2 . The three 

couplings in (3.12) are related to the single parameter g 2 : 

9 D U D U* = - 2 92^Hm DU 

Od^d^ = 2g 2 ^m^m D(s) m D * s) (3.13) 

9d (s) d (s) ^ = 2g 2y /m^m D{s) . 

In principle, the couplings g\ and g 2 must be computed by nonperturbative methods. An 
estimate can be obtained invoking vector meson dominance (VMD) arguments. For example, 
one can consider the D-meson matrix element of the scalar cc current: (D(v')\cc\D(v)) , 
assuming the dominance in the t-channel of the nearest resonance, i.e. the scalar cc state, 
and using the normalization of the Isgur-Wise form factor at the zero-recoil point v = v'. 
This allows to express Qddxco m terms of the constant f Xc0 that parameterizes the matrix 
element 

(0\cc\xM) = /xco™xco , (3-14) 

obtaining: 

9dd X co - 1 — 7 > [6. Lb) 

JXcO 

a relation which determines g\ once f Xc0 is known: 
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Adopting the same argument one can also obtain g 2 in terms of the J/ip leptonic constant 
fy, defined by (0|c7 M c| J/ip(p, e)) = f^m^. From the VMD result 



9DDi> = ~r (3.17) 

u 



one gets: 



92 = • (3.18) 



The input quantities for computing the diagrams in fig.l are now available. We have 
only to notice that the strong couplings described above do not account for the off-shell 
effect of the t-channel Z^*j particles, the virtuality of which can be large. As discussed in 
the literature, a method to account for such effect relies on the introduction of form factors: 

9i(t)=9 i0 F i (t), (3.19) 

with g i0 the corresponding on-shell couplings (2.9), (3.7) and (3.9). A simple pole represen- 
tation for Fi(t): 

F t (t) = A ^ gf } (3.20) 

is consistent with QCD counting rules [29]. We adopt it, keeping in mind that the parameters 
in such form factors represent a source of uncertainty in our results. 



IV. NUMERICAL ANALYSIS OF B~ -► K~h, 



c 



Considering the diagrams in fig.l with M Sc = h c , there are ten possible combinations of 
intermediate states corresponding to non-vanishing strong vertices. Some of such diagrams 
vanish, since the rescattering amplitude is parity conserving and the final state K~h c has 
positive parity due to angular momentum conservation. As a consequence, only the parity- 
violating weak decay amplitude contributes, hence only the intermediate states (D s , D) and 
(D* a ,D*) must be considered in fig.l. The expression of the absorptive part of a generic 
diagram reads as: 

,/A(m 2 B , m' 2 ( , b mL,,) r +i 
imA = V ^ D()) f d zA{B- -> D^-D^°)A(D^-D^° -> K~h c ) . (4.1) 

In the case of the diagram in fig.l corresponding to B~ — > D~D° — > K~h c mediated by D*° 
the expression (4.1) becomes: 

_ K^/m B m D o v 
im^li - — 327rm 2 — A {m B ,m D -,m D o)j Ds 

V 2rriBm D o I J-i t — m z D , 
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with K = ^V c bV* s ai, A the triangular function, q the kaon momentum and e the h c polar- 
ization vector. The function f 1 is given by: 



AC*) = - 



m D J m D 



{( ml - q - k -l 

V m D* 



with q o = nBTH'K-H'h, - = 2_ ^B,"^,"^ ^0 

2tob 2rriB 2uib 

/\ I 771 d ; 771 r-jo ; 777- n J 7 * f) i — *i i i 9 9 * th 

— , q • k — q k — \q\ \k\z and t = m K + m D — 2q ■ k. Expressions 

for the other diagrams can be worked out, analogously. The t-dependence of the couplings 
is given by eq.(3.20) with all Aj put equal to a unique parameter A. 

We use \V cb \ = 0.042 and |V CS | = 0.974, the central values reported by the Particle Data 
Group [22], and ai = 1.0 as obtained from the analysis of exclusive B — > £)(*)£)(*) transitions 
[25]. Exploiting the heavy quark limit, we put fn* = /d s and use fo 3 = 240 MeV [7]. As 

for the Isgur-Wise form factor, the expression = (j^j is compatible with the current 

results from the semileptonic B — > decays, and the product V c ^ coincides with the 
experimental determination reported in [25]. 

A comment is in order about the vertices D^D^K, expressed in terms of the coupling 
g according to (2.10). An experimental determination of g has been obtained by CLEO 
Collaboration measuring the full D* width and the D* branching fraction to Dn. The 
result is g — 0.59 ±0.01 ±0.07 [30]. Such a determination should be compared to theoretical 
predictions ranging from g ~ 0.3 up to g ~ 0.77 [31]. Since the expressions of the rescattering 
amplitudes always contain the product of g and the form factor (3.20), we use the central 
value of g obtained by experiment, leaving to the parameter A the task of spanning the 
range of possible variation of the coupling. 

For g 1 we use eq. (3.16) together the QCD sum rule result f Xc0 = 510 ± 40 MeV [12]. 
The coupling #2 can be obtained using (3.18) and the experimental value /;u = 405 ± 14 
MeV. The VMD determination of the J/ip couplings is reproduced by QCD sum rule and 
constituent quark model analyses [32]. Relating the various couplings to g 1 and g 2 we use 
m D = m Ds and m D * = m D *. 

Eq. (4.1) allows us to compute the imaginary part of the rescattering diagrams. 
The determination of the real part is more uncertain. A dispersive integral may be 

used: ReA(m R ) = — PV / — — l —^-ds with the thresholds s+J given by: s\ l J = 



^ 2 Q ' k - m 2 K )k° - (m B - q°)q ■ k}} (4.3) 

m z D , m B \q\ z 

+ m 2 K -m 2 hc X^ 2 (m 2 B ,m 2 K ,m 2 hc ) m\ + m 2 Ds - m 2 D0 , r . 



[m 



D 



(*) +m D (,)) 2 for any specific diagram. Assuming that the integrals are dominated by the 

s 

region close to the pole m 2 B , so that they can be computed by using a cutoff not far from the 
B meson mass, we obtained for B~ — > K~Xco{J/ip) that the real parts of the amplitudes 
are approximately equal to the imaginary parts, with large uncertainties due to the cut-off 
procedure [12]. For this reason we account for the real part of the amplitudes considering 
them as fractions of the imaginary part varying from to 100%, i.e. we include their con- 
tribution to the final result considering the range from Re^4j = up to Re^4j ~ InuAj. Such 
an uncertainty cannot be removed in our approach and will affect the final result. 
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A parameter is left in our analysis, i.e. the constant A in the form factors (3.20). One 
would expect A of the order of the mass of radial excitations of the charmed mesons. It 
is possible to constrain the range of values for such a parameter considering rescattering 
contributions to B~ — > K~J/tp, where the sum A(B~ — > K~J/ip) = A fact + A resc is 
bounded by the experimental measurement of the branching fraction B(B~ — > K~J/ip). If 
one considers the range 2.6 — 3 GeV for A one gets a rescattering contribution not exceeding 
the experimental bound. Moreover, one can consider B~ — > K~XcO as provided only by 
rescattering effects, repeating the analysis in [12], with the difference of using the relations 
(3.9) which imply a factor of 3 between the couplings of Xco to pairs of D and D* mesons, 
dictated by the spin symmetry. With this factor into account, one gets a branching fraction 
compatible with the experimental result from BABAR if the parameter A is varied around 
3.0 GeV. 

Provided with such constraints we analyze B~ — * K~h c . In fig. 2 we plot the branching 
ratio obtained considering the rescattering amplitudes as a function of A. We find a region 
that can be represented by the interval: 

B(B~ -> K~h c ) = (2 - 12) x 10~ 4 , (4.4) 

where the range of values accounts for the uncertainty on the dispersive part of the rescat- 
tering amplitudes and on the variation of the parameter A. This result suggests that 
B~ — > K~h c occurs with a rate large enough to produce a signal at the B-factories, as 
discussed in the next Section. Moreover, the outcome (4.4) implies that B~ — > K~h c 
represents a sizeable fraction of the inclusive B~ — > Xh c mode, the branching ratio of 
which, estimated considering the production of the cc pair in h c in the color-octet state, is: 
B(B~ -> h c X) = (13 - 34) x 10- 4 [33]. 

The theoretical uncertainties affecting our results are related to the poorly known values 
of some of the input parameters and to the basic assumptions adopted in the calculation. 
While the numerical values of several parameters (namely, the strong couplings among heavy 
mesons) can be made more precise using new experimental or theoretical information, it is 
difficult to assess the actual size of the uncertainties related to the computational scheme we 
have used in evaluating rescattering effects. The main uncertainty in the numerical results 
is due to large cancellations between different amplitudes, which individually turn out to be 
of similar size. This is common to calculations involving hadronic degrees of freedom, and 
it is not easy to envisage a procedure for reducing or controlling the final error. Another 
uncertainty is due to the neglect, in the calculation of diagrams in fig.l, of contributions of 
higher resonances and of many-particle intermediate states, even though a minor role can be 
presumed for higher resonances since the corresponding universal form factors and leptonic 
decay constants are expected to be smaller than for low-lying states. 

Bearing such uncertainties in mind we can conclude that rescattering terms may con- 
tribute to the nonfactorizable effects observed in B — > charmonium transitions. 

V. REMARKS ABOUT THE OBSERVATION OF B~ -► K~h c AND 

CONCLUSIONS 

Let us discuss few phenomenological consequences of our study, coming first to the 
possibility of detecting and studying h c using B decays. 
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As mentioned in the Introduction, observation of h c has been reported in pp annihilation 
and in p — Li interactions, where the meson is produced through qq annihilation in three 
gluons. Other production mechanisms are possibile at e + — e _ machines, namely via ip' 
intermediate production. For example, one can consider the radiative decay ip' — > 77^,7 with 
the subsequent transition i]' c — > h c ^ as feasible to obtain a sample of h c . Another possibility 
is the hadronic decay mode ip' — > h c 7r°. In this case the estimated branching ratio is rather 
sizeable: B(ip' — > h c 7i°) ~ O(10~ 3 ) [34], and therefore one could consider the investigation 
affordable, e.g., at a charm factory; however, a low n° reconstruction efficiency could severely 
limit the possibility of studying h c produced by this decay chain. 

As for h c produced in B decays, one could access the meson looking either at its hadronic 
modes: h c — > J/tpn°, p°7r°, /ii/ (980), h\KK, . . . , or at its radiative modes: h c — > ^7, Xcol, 
etc. In particular, the channel h c — > rjc'j seems promising, as noticed by Suzuki [35]. Its 
branching ratio, estimated assuming that the h c wave function close to the origin is the same 
as that of Xcu is large: B(h c — > 77^,7) ~ 0.50±0.11 [35]. A similar result: B(h c — > rjc'j) = 0.377 
[36] is obtained using the charmonium wave functions parameterized in ref. [37]. These two 
predictions, together with the experimental datum for B{r\ c — > KKir), allow us to translate 
our result (4.4) in a prediction for the decay chain B~ — > K~h c — > K~r\ c ^ — > K~ (KKn)j 
which can be studied at a B factory: 

B(B~ -> K~h c -> ir-r/ c7 -> K'iKKn)^) = (4 - 26) x 10~ 6 , (5.1) 

a result within the reach of current experiments. It is worth noticing that the investiga- 
tion of this particular decay chain is favoured by the rather accurate knowledge of the r\ c 
hadronic decays, and by the fact that one could use the r) c mass and the photon direction 
to discriminate the signal from the background. 

Coming to the role of rescattering effects in B — > charmonium transitions, we have found 
that they can be effective, and are able to produce for the mode B~ — > K~h c a branching 
fraction comparable with that of B~ — > K~Xco- Further evidence for the presence of large 
nonfactorizable contributions in B decays with charmonium in the final state can be obtained 
by looking at other decay modes. One possibility is B~ — > K~ip(3770) which, because of 
the smallness of the leptonic decay constant ^,(3770) , is predicted by the factorization model 
with a tiny branching ratio. The observation of this decay mode with a sizeble branching 
fraction B(B~ -> 7^(3770)) = (0.48±0.11±0.12) x 1(T 3 [38] represents a further evidence 
of the presence of large nonfactorizable contributions. In our approach, using g_DD</>(3770) — 
14.94 ± 0.86 obtained from the width of ^(3770), we would get B(B~ -> ^"^(3770)) = 
(0.9 — 4) x 10~ 4 , consistent with the experimental datum considering the large theoretical 
uncertainty. Similar conclusion applies to B~ — > K~Xc2 with Xc2 the J PC = 2 ++ state of 
the charmonium system, the amplitude of which also vanishes in the factorization approach. 
The observation of this decay mode with branching fraction comparable to B(B~ — > K~Xco) 
and B(B~ — > K~h c ) would support the rescattering picture. 
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FIGURES 




FIG. 1. Typical rescattering diagrams contributing to the decay B~ — » K~M c5 , with M C c a 
meson belonging to the charmonium system. The boxes represent weak vertices, the dots strong 
couplings. 




FIG. 2. Branching fraction B(B — » h c ) versus the parameter A. The lowest curve corre- 
sponds to ReAi = 0, the highest one to ReAi = ImAi. The dark region corresponds to the result 
(4.4). 
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